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Abstract 

The inclusive semileptonic branching ratio is often employed to normalize other 
inclusive B decays. Using recent determinations of the non-perturbative parameters of 
the Operator Product Expansion we compute the normalization factor for the branch- 
ing ratio of B — > X s j and find a few percent enhancement with respect to previous 
determinations . 



1 Introduction 



The partial width of the inclusive B decays to light quarks depends on the fifth power of the 
b quark mass. In order to reduce the uncertainty that follows from this strong sensitivity 
and to avoid the large radiative corrections that are sometimes related to the heavy quark 
mass, the Branching Ratio (BR) of rare decays is usually expressed in terms of the CKM 
favored semileptonic BR, BR C &, = BR[£> — > X c £u], a quantity that is presently known at 
the 1% level [U [2]. This choice introduces a marked dependence on the charm quark mass 
in the calculation of rare decays, but it is very convenient in b — > s transitions, whose 
CKM factor |Vf*V^| is essentially determined by \V cb \ measured in inclusive semileptonic 
B decays. Moreover, in the case of B — > X^ the charm mass dependence due to the 
normalization to BR C ^ is partially compensated by that of the perturbative QCD corrections. 
In recent years the B factories have performed increasingly detailed studies of semileptonic 
B decays, providing us with improved determinations of the b and c quark masses and of 
the Operator Product Expansion (OPE) parameters. In parallel, both the measurements 
and the theoretical calculations of inclusive rare decays have improved significantly. In this 
Letter we reconsider the normalization of rare decays and try to assess its uncertainty, taking 
into account the latest developments. We will concentrate on the radiative inclusive decay of 
the B meson, B —> X s '~f: the Next-to-Next-to-Leading Order (NNLO) calculation of its BR 
in the Standard Model is quite advanced [3] and its experimental error will soon approach 
5%. Many of our considerations apply to B — > X s £ + £~ j3] and B — > X u £v as well. 

As mentioned already, the unitarity of the CKM matrix implies |V^*V^| 2 = [1 + A 2 (2p — 
l)+0(X 4 )}\V cb \ 2 = (0.963±0.003) \V cb \ 2 [5], i.e. the CKM prefactor of B -> X s7 is essentially 
given by \ V c b\, whose best determination follows from BR^,, and has a 2% error pp. In recent 
analyses [H [7] the radiative BR is therefore expressed by 

BR 7 (^ ) = dK[B -> X sl \ E >Eo - I — 2 -— 

where the phase space ratio C is defined by 
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T[B^X u ev\ (2) 

In Eq.([T]) the radiative width is normalized to the charmless semileptonic decay in order to 
split the charm mass dependence of the perturbative matrix elements of b — > X s 7 — an 
0(a s ) two-loop effect — from that due to the normalization — a tree- level effect, contained 
in C |6j. This choice makes the calculation more transparent. One has 



T[B -> X sl } Ei>Eo 
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where a em is the fine structure constant. The perturbative term P(E ) is estimated at NNLO 
in QCD P, [7] and includes electroweak effects [8], while Snp contains the non-perturbative 
power corrections [9]. Fortunately, the m c dependence of P(E ), despite being a loop effect, 
compensates about half of that of C. 
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2 Calculation of C 



Like all inclusive widths, the ratio C can be calculated using the OPE and expressed as a 
double expansion in a s and inverse powers of the b quark mass, currently known through 
0(a 2 ) and 0(Aq CD /ml). C depends sensitively on the b and c quark masses, as well as on the 
matrix elements of the dimension 5 and 6 operators. This is where the recent experimental 
studies of the inclusive moments of B — > X c ev and B — > X s 7 enter in a crucial way. 
Indeed, the moments of various kinematic distributions provide information on the non- 
perturbative parameters of the OPE. Global fits to the moments describe successfully a 
variety of moments and allow for a 40 — 50MeV determination of m c and mj, a ~ 10 — 20% 



w, 



cb\ 



determination of the 1/m 2 and 1/mf matrix elements, and a ~ 2% determination of 
[21 [TO]. There are different ways to take into account the available information, relying on 
different assumptions and schemes. We work in the kinetic scheme where a 'hard' cutoff 
H separates perturbative and non-perturbative effects respecting heavy quark relations, and 
non-perturbative parameters are well-defined and perturbatively stable. 

Our starting point are the NNLO expressions for the charmed and charmless total 
semileptonic widths 
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where a s = ai ri/_5 ' ) (mb), r = (m c /m;,) 2 , g{r) = 1 — 8r + 8r 3 — r 4 — 12r 2 lnr, and all the masses 
and OPE parameters are defined in the kinetic scheme at finite m& with /i ~ lGeV. The 
non-perturbative corrections have been computed in [12] and are expressed in terms of the 
parameters /i 2 . 
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p\ s . The matrix element of the Weak Annihilation (WA) operator 



B WA = (5|0^ A |-B) is poorly known. It is here renormalized in the MS scheme at the scale 
/i WA , see [321 E]. We recall that B WA vanishes in the factorization approximation, and that 
WA is phenomenologically important only to the extent factorization is actually violated. 
There is however an 0(1) mixing between WA and Darwin operators, and at lowest order 
in perturbation theory one has _B WA (//) = B WA (fi) — p 3 D /2n 2 Inp'/fi. As factorization may 
hold only for a certain value /i WA = [if for which B WA (fif) = 0, a change of the scale /i/ 
provides a rough measure of the (minimal) violation of factorization induced perturbatively. 
We neglect intrinsic charm contributions [15] . WA uncertainties make a precise prediction 
of C problematic at present. Fortunately, they cancel out in Eq.Q since the radiative BR 
cannot depend on the non-perturbative features of the charmless semileptonic decay. 
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The perturbative corrections at \i 



(on-shell scheme) are given by 
25 2 „ 
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where h(r) [16J is given in App. C of |6j, and (3 { 
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11 — |n/. The O(a^) perturbative 



corrections are known exactly in both cases [XT] 18J. Their m c dependence is given here 
in terms of simple interpolation formulas, valid for m c /m b between 0.2 and 0.3. While the 
0{a 2 s (3o) part of p^ has been known for some time [TH], the remaining (non-BLM) term 
is a very recent result [18J. The /i-dependence of Pu,' c (r,p,) can be found exploiting the 
//-independence of the widths at each perturbative order and the known /i-dependence of 
masses and OPE parameters [20] 
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In the above formulas Cf = §, CU = A c = 3, and we have used n/ = 3 in /3 , corresponding 
to three light massless quarks: indeed, the m c dependence of the charm loops on gluon lines 
is not known, and one can approximately decouple them. This is also consistent with the 
calculation of the moments in the kinetic scheme [21]. A fully consistent 0{a1) implemen- 
tation of the kinetic scheme would require the 0(a s ) corrections to the Wilson coefficients 
of the higher dimensional operators in (4|5), which is not yet available except for the trivial 
/4 term. As for higher order power corrections, we recall that 1/77t| corrections have been 
estimated to be tiny in the charmed decay rate [22j . 

The numerical value of C depends on those of the OPE parameters and of a s . We take 
as default values the results of the global fit [2], namely 
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where all values are in the kinetic scheme with p, = lGeV. We employ a s (m b ) 
set /i WA = m b /2, obtaining 



0.219 and 



C = 0.546 - 2.0 B WA (m b /2) 
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In the second line we have listed the individual contributions, that are accidentally all 
negative. Using only BLM corrections and setting n/ consistently equal to 3, one gets 
C = 0.543 — 2.0B WA (rrib/2), which shows that the 0(a 2 ) corrections are dominated by 
O(a 2 j3o) running coupling effects. The 0{a 2 j3o) corrections can be absorbed in a rescaling of 
a s in the NLO contribution, using ot s (pblm) • It is clear however that in that case the BLM 
scale is very low, just above 1 GeV. Following [U] we use B WA (m b /2) ~ as central value 
and vary it between and +0.012 GeV 3 , since the positivity of the q 2 spectrum of B — > X u tv 
at p WA = nib/ 2 suggests a positive WA contribution. Taking into account the correlations 
among the input parameters [2], and varying /i between 0.7 and 1.3 GeV to estimate the 
residual perturbative error, we find 

C = 0.546 ± 0.017(par) ± 0.016(pert)+°;°°°(WA) = 0.546^|. (7) 

This value can be compared with C = 0.580 + 1.8 B WA (nib) ±0.016, obtained in a global fit 
to semileptonic and radiative moments [TU] where the b mass is in the 15 scheme [23] • Using 
the value of p 3 D (fi = 0) obtained in [TP] , we can rewrite it as C = 0.582 + 1.8 B WA (nib/2), that 
differs from Eq.([6]) by 6.5%, exceeding the stated uncertainties. We observe that the fit in 
[TO] is based on older data than that of [2]. The two fits also differ in the perturbative scheme, 
in several assumptions, and in the estimate of theory errors, as detailed in Refs. [10] and 
[2"T] 12"]. resp. Ref. [TU], for instance, does not extract the charm mass directly but eliminates 
it through the heavy meson mass relation 

m b ~ m c = M B - M D + & ( — - — ) + (J--l-\ + 0(l/m 3 ), (8) 

2 \m c m b J 4 \mi m^ J v 

where all quark masses and OPE parameters have the same normalization point p. This 
relation involves an expansion in l/m c , rather than 1/nib, and some non-local operators 
(the term p 3 ) that do not enter the expressions for the semileptonic widths. In [TUJ the 
normalization scale is set to zero and the fit is effectively equivalent to a direct fit to the pole 
mass difference nib — "V The default procedure in [lU] employs other meson mass relations 
to fix \Xq and p| 5 , treats differently A terms, includes the moments of B — > X s j accounting 
for neither distribution function effects nor an additional theory error, and estimates theory 
errors in a different way. 

Despite their differences, the two fits give compatible values of | V^t, | . But the ratio C, well 
approximated by 1.2 — 2.2 m c /nib in the physical region, is more sensitive to the exact value 
of the charm and bottom masses than \ V c b\, due to the correlation between nib and m c in the 
results of the fits. As illustrated in Fig. 1, the semileptonic rate and moments (in particular, 
those of the lepton energy distribution) stringently constrain a certain combination of m Ct b, 
given by m nib— 0.65 m c that roughly corresponds to constant values of the semileptonic total 
width for fixed \V c b\. In the fit of [2] the latter gets a 0.3% relative parametric error, while 
C ~ \.2 — 2.2m c jnib has a 3.3% error. The charmed inclusive width, roughly proportional to 
ml (1.2 — 2.2m c /mb), has instead a 1.2% relative uncertainty. As shown in Fig. 1, while the 
fitted value of ml (1.2 — 2.2 m c /m b ), and consequently of \V c b\, is insensitive to small changes 
in the data and to the inclusion of radiative moments, the situation is quite different for m c 
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Figure 1: la contours in the m c , m b plane in the kinetic scheme at ix — 1 GeV from: a) the global 
fit |2J (small red ellipse in the center); b) same fit without the radiative moments (larger upper 
right ellipse); c) the Belle fit |24| (larger lower left ellipse); d) the PDG [25] (large, light green 
central ellipse); e) the sum rules determinations [26J (smaller, light red, center bottom ellipse), 
after translating them in the kinetic scheme. The solid red and dashed blue lines correspond to 
constant values of the semileptonic width and of C , respectively. C decreases moving to the 
right of the plot. 



and C. As indicated by the IS* fit to only Belle data performed by the Belle collaboration 
[24£j the treatment of theoretical errors may similarly have a larger impact on the quark 
masses than on \V c b\. Different determinations of the c and b masses tend to prefer values 
close to the center of the plot: we show in Fig. 1 the la regions in the (m c , mb) plane selected 
by the PDG [25] and by the a(e + e~ —>■ hadrons) sum rules according to [26] , taking into 
account the non-negligible error introduced by the scheme translation (40 and 50 MeV for 
mb and m c , respectively, estimated using the residual scale dependence). 

It is instructive to compare results based on the same set of data, like the recent Belle 
fits [24J that follow the methods of [5TJ and [TU]. Using their kinetic scheme results, in 
the same way that has led to Eq.([6]), we obtain C = 0.574, although with a parametric 

1 The Belle data are also included in the global fit of [2] • 
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error twice larger than above. From the central values of the IS scheme fit of [24J we find 
ml° le — m p c ole = 3.393 GeV that can be employed to eliminate m c in Eq.Q, using the IS 
scheme for m&. This is analogous to what was done in Appendix C of [6], but includes 
l/m| effects that were neglected in that paper. Assuming B WA (rrib/2) = 0, the result is 
C = 0.563^] In other words, using the same data set the results in the two schemes are 
closer and the kinetic scheme gives a higher value. The discrepancy between Eq.(|6]) and 
[TU] is likely to be due to different input data, although at some level differences in the fit 
methods and in the treatment of theory errors may also play a role. The situation is likely 
to improve with better experimental data and better control of theory errors (an upgrade of 
the fitting methods is under way [27]). but in the meantime the error attached to C should 
be treated with caution. 



3 Impact on radiative decays 

We now move to the calculation of the radiative BR. The power correction Snp to the total 
inclusive radiative decay normalized to B — > X u £u [9] — see Eq.pl) — is 



>NP 



\27rn;? 21ra 2 c m h ) Cj°^ efl (fib) V ^ m 1 J m l m l 



where /i^ is the low-energy scale, that we set equal to 2.5 GeV like in [7J, C^^ifi 
—0.37 is the effective Wilson coefficients of the operator Qj calculated at leading logarithmic 
order, and C^ \fib) ~ 1.20 is the appropriate combination of Wilson coefficients of Qi >2 
(see P I7J)- An additional 0(a s /i^/m^) contribution, dependent on the photon energy cut, 
has been calculated in [28], but it is small for E < 1.6 GeV and we neglect it. In the 
notation of [7J, 8 NP = N (E ) / P (0) (ji b ) . Since the NNLO calculation of P(E ) P Cj has 
not been performed in the kinetic scheme, all the parameters in Eq.([9]) must be understood 
at ii — 0, but we convert them to /x = 1 GeV expanding to second order in a s . We find 
$np — 0.033 — 3.2B WA (rrib/2). Since 6^p is correlated with C, a quantity that summarizes 
all power corrections and the normalization in Eq. (1) is 



b 



1 + 5 



j pIs 



= 1.600 + 0.075 Qs + 0.061 Q 2 + 0.047^ + 0.057^ + 0.019, 
= 1.859 ± 0.054(par) ±0.045(pert), (10) 

which is independent of both B WA and /i WA . As long as the perturbative corrections to Eq.([9| 
are not available, the numerical value of the charm mass to be employed is quite arbitrary. 
We have employed our default input value, but using m c ~ 1.5 GeV, closer to the pole mass, 



would decrease F by 1%. We have accordingly increased the perturbative error in Eq.(10). 
We list in the Appendix approximate formulae that allow for an easy determination of C and 
F with different inputs. In the presence of new physics C^ eS (/j,b) = C^gM (fib) + ^7new (/■*&) 



^Employing the same method with the results of [TU], we get C = 0.585. 
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and this is expected to be the only change in F, leading to 



F = F SM - (0.06 ± 0.03) e new (11) 



with e new = C^J(// 6 )/Cj 0)eff (/i 6 ) and F SM given in Eq.(10). 

The BR with a photon energy cut at the standard value E = 1.6 GeV can be computed 
from Eq.Q using the approximate relation 



P(1.6 GeV) ~ 0.1247 - 0.0572 [m c (m c ) - 1.224 GeV] + 0.0084 {m\ b - 4.68 GeV) (12) 

that can be extracted from [7J. Since this result is obtained with in the IS scheme 
and m c in the MS scheme, we perform an explicit change of scheme to the kinetic one, 
differentiating the NLO contribution to P(Eq) wrt m c ^. Inserting Eqs.(ToJT2) into Eqs.(3) 



and (1) and using BR C ^ = 0.1064 from [2J, one gets BR 7 = 3.28 x 10 4 which we choose as 
our central value. Alternatively, we can convert our kinetic scheme mass inputs to the IS and 



MS schemes and employ the results in Eq. (12), facing however an extra theoretical error due 



to the conversion. As both masses enter P(E ) at 0(a s ), one can use one-loop formulas with 
a s (2.5GeV) « 0.27, namely use m\ s = 4.69 GeV and m c (m c ) = 1.23 GeV in EqJT^). This 
leads to BR 7 = 3.30 x 10 -4 . If the conversion of the kinetic masses to the IS and MS schemes 
is made via two-loop expressions, BR 7 can be as low as 3.21 x 10~ 4 , depending on the scale 
chosen for the m c conversion. The spread of these values is an estimate of the perturbative 
higher orders consistent with that given in [3]. All renormalization scales involved in the 
calculation of P{Eq) are kept to the default value of [7j. The parametric error due to the 
OPE parameters and BR C ^ given by the fit follows from the correlation matrix in [2] and 
amounts to only 1.8%, to which we add in quadrature the uncertainty from the CKM factor, 



a s , m t , and the theory error in Eq.(10), obtaining 3.7%. Summarizing, we have 



BR 7 (1.6 GeV) = 3.28 x 10~ 4 [1 ± 0.037 ± 0.03 ± 0.03 ± 0.05] , (13) 

where the four errors are due to i) the normalization and parametric uncertainties; ii) the 
perturbative uncertainty in P(1.6GeV); in) the m c interpolation of [7j; iv) unknown non- 
perturbative contributions beyond the OPE. Apart from i), we have employed the same 
errors as in [31 [7] and our central value is about 4% higher than there. The main reason 
for the shift is the new determination of C . Refs.[3l [7] employed the most precise value of 
C available at that time [TU]. Our total 3.7% normalization and parametric error is a little 
larger than the corresponding uncertainty given in [3], [7] , but the total error is still dominated 
by the ±5% non-perturbative uncertainty. 

The slow convergence of the perturbative series in Eqs.(6jl0) is slightly disturbing. It is 



partly accidental, as both numerator and denominator have ~ 2% second order perturbative 
corrections, but with different sign. On the other hand, the kinetic mass definition is not 
particularly appropriate for a quark with mass m ~ fi. An alternative is to use a hybrid 
scheme where m;, and the non-perturbative matrix elements are defined in the kinetic scheme, 
while the charm mass is defined in the MS scheme. This yields a better apparent convergence 
in both C and T[B — > X c ev] if one employs m c (m c ), namely the MS scale of the charm 
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mass is set equal to the mass itself. As for the numerical value of m c (m c ), to simplify the 
comparison with Eqs.( 6fTo" ) we choose m c (m c ) = 1.267 GeV that is obtained from our input 
in (Jsl) through the two-loop perturbative relation with a s (mb)- The results are 
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0.574 - 0.005 Q 
1.886 



,(m 6 /2) 



0.001 a 2 - 0.004^ - 0.022^ - 0.001 pis - 1.9 5 WA ( 



m b /2) 



(14) 



1.741 + 0.019 Q , + 0.003 a2 + 0.047 



0.059„3 +0.018„ 3 



which agree with but converge better than Eqs. fl6fl0] ). Since direct fits to m c (m c ) and 
kinetic scheme OPE parameters will soon become available [27], we list in the Appendix 



approximate formulae for this option too. The /^-dependence in Eqs. (14) is less than 0.4%. 



One can also calculate C and F using the 15* definition of the bottom mass [231 and the 
MS scheme for m c : using again inputs that correspond to those employed in Eqs.( 6p0 ), i.e. 



AS 



= 4.74 GeV and m c (m c ) = 1.267 GeV, we find C = 0.545 and F = 1.870. Our results 
are quite stable for changes of scheme. 

Finally, let us consider alternatives to the normalization of Eq.(l). Since much of the m c 
sensitivity is related to the normalization to BR C ^, one could hope to reduce it by avoiding 
BR C &,. This amounts to reexpressing the first ratio in Eq.(l) in terms of the semileptonic 
charmless width: 
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where tb = 1.585(7) ps is the lifetime of the B°/B ± admixture and T[B — > X u ev\ is given 
in Eq.([5]). As \V u b\ drops out, the input parameters necessary to compute the rhs of Eq.(15) 
are |V^>|, m^, and the OPE parameters. All these parameters, as well as m c and BR C &,, 
can be extracted from a global fit as in [2J. Clearly in that case the radiative BR will be 
the same and, because of the various correlations, will have the same error if computed 



using the rhs or the lhs of Eq.(15). The rhs is very sensitive to due to the m\ factor in 



Eq.([5]) but insensitive to m c (which however still enters P(E )). Moreover, the theoretical 
uncertainty in the extraction of \V c b\ from BR^ affects only the rhs of (Jl5|J^l Neglecting all 
the correlations between the OPE parameters and BR^ 
5% parametric error instead of 1.8% 



get a 



and using the rhs of Eq.(15), we 
Even the very small errors on m c ^ found in 
would lead to a 3.3% uncertainty in this case. This demonstrates the advantage of the 
normalization to the semileptonic BR and of a global fit to the OPE parameters, as already 
stressed in [7J. In the future, independent and precise determinations of and m c could 
be used as additional constraints in the fit. 



4 Summary 

We have calculated the normalization factor for radiative inclusive B decays in the kinetic 
scheme and discussed its uncertainty and dependence on the input parameters, taking into 

3 See eq. (7) in the second paper of Ref. [2D] for an estimate. 



S 





1 1 tin DljllClllC 








C G 




C LS 




c 


kin 


0.9185 


0.035 


-0.001 


-0.021 


-0.186 


0.005 


-0.53 


F 


kin 


1.085 


-0.045 


0.007 


0.148 


0.169 


-0.091 


0.60 


C 


MS 


1.001 


0.029 


-0.102 


-0.021 


-0.186 


0.005 


0.03 


F 


MS 


1.001 


-0.035 


0.112 


0.148 


0.169 


-0.091 


-0.01 



Table 1: Coefficients of the approximate formulae for C and F for different m c schemes. 



account the correlations that arise from the fit to the moments of semileptonic and radiative 
distributions. Using the latest global fit in the kinetic scheme [2] we obtain C = 0.546+o;q33. 
We estimate that adopting our normalization in the first NNLO analysis of [31 Ej would lead 
to a 4% higher BR 7 for £ 7 > 1.6 GeV, (3.28 ± 0.25) x 10" 4 . Since the normalization factor 
is more sensitive to the value of the c and b quark masses than \V c b\, its determination is 
presently more volatile, but progress will come from improved measurements of the moments, 
a better understanding of the theoretical uncertainties involved in the fits, and complemen- 
tary constraints on m C) &. First steps to reduce the theory error in the normalization would 
be to employ the MS charm mass and to compute higher order perturbative corrections, 
starting with the implementation of [18] . 
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Appendix 



The factors C and F defined in Eqs.(2jl0) are well approximated by the following formulae, 



whose coefficients are given in Table 1. 

327T 2 

ci + c b 5 b + c c 5 C + c G [m g + c D p D + c LS pis + c as 5 as - -^B(m b /2) 



C = g{r 

F = \g(. r )Y l [ci + c b 5 b + c c 5 c + c D p 3 D + {c G p G + c LS p\s) /m 2 c + c a J aa ] 



m l 



Here all parameters are in the kinetic scheme with p = 1 GeV, except for the charm mass that 
is either in the kinetic scheme or in the MS scheme, m c (m c ). Moreover, 5 b = m b — 4.6 GeV, 
6 C = m c - 1.15 GeV, and S aa = a s - 0.22. Notice that in the range 0.22 < y/f < 0.29, 
g(r) ~ 1.1928 — 2.2443 m c /m b within 0.03%. The approximate formulae have a precision 
better than 0.3% in the ranges 4.5 < m b < 4.7 GeV, 1 < m c < 1.3 GeV, 0.2 < a s (m b ) < 0.24. 
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